Signature of QCD critical point: Anomalous transverse velocity dependence of 

antiproton-proton ratio 
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We formulate the QCD critical point focusing effect on transverse velocity (fit) dependence of 
antiproton to proton (p/p) ratio, which was recently proposed by Asakawa et al. as an experimental 
signature of QCD critical point in high energy heavy ion collisions (HICs). For quantitative analysis, 
Ultra-relativistic Quantum Molecular Dynamics (UrQMD) transport model and THERMal heavy- 
IoN generATOR (THERMINATOR) are applied to calculate the corresponding j3 t dependence of 
p/p ratio for three gedanken focused isentropic trajectories with different focusing degree on QCD 
phase diagram. Finally, we obtained an observable anomaly in /3 t dependence of p/p ratio, which 
can be employed as a signature of QCD critical point. 
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I. INTRODUCTION 



In recent years, ultra-relativistic high energy heavy 
ion collisions (HICs) experiments, such as SPS/CERN 
(v/snn ~ 10 A GeV) and RHIC/BNL (^/s^ ~ 200 A 
GeV), have been aimed to search for the new form of 
matter, which is composed of deconfinement free quarks 
and gluons, and thus called Quark Gluon Plasma (QGP) 
[2 US El- Lattice-QCD calculation predicts that the 
phase transition between hadronic and QGP phases at 
vanishing baryon chemical potential ps is crossover tran- 
sition whereas at higher fiB the transition may become 
first-order 0, Q. However, there are many uncertain- 
ties for Lattice-QCD or model calculation to determine 
the first-order phase transition boundary as well as the 
location of corresponding end point, the so-called QCD 
Critical Point (QCP) 0. Thus, the location and even the 
existence of the QCP are still open questions. The uncer- 
tainty of theoretical calculation drives us to explore the 
properties of hot QCD matter at higher net baryon den- 
sity and find experimental evidence of the QCP. Recently, 
several experimental programs have planed to search for 
the QCP, such as the energy scan program of RHIC/BNL 
(V s nn = 5 ~ 39 A GeV) [a, EJUS], the light ion program 
of NA61 experiment at SPS/CERN (y^f = 5 ~ 17.3,4 
GeV) [HEa] and also the CBM experiment at FAIR/GSI 
(a/snn < 8.5A GeV) [H, 0, EH . Many experimental ob- 
servables have been also proposed to be the QCP signa- 
tures, such as dynamical fluctuations in K/tt ratio [16j |. 
two experimental observables correlation [ljj , high order 
moment of transverse momentum [l8j . etc.. 

In the hydrodynamical description of relativistic heavy 
ion collisions, the expansion of the central fireball should 
be regarded as isentropic for the negligible entropy pro- 
duction in the latter evolution, which can be represented 
as entropy density s to baryon density rif, ratio (s/rib) to 
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be constant [191 ]. A trajectory with s/rif, = constant for 
a given colliding system on the QCD phase diagram is 
called isentropic trajectory. Recently, the critical singu- 
lar properties of QCP have been implemented in equation 
of state (EOS) for hydrodynamical description of HICs 
[ISIj [2QJ by Asakawa, et at. They pointed out that when 
the isentropic trajectory passes through the vicinity of 
the QCP, it may be deformed-the so-called "QCP focus- 
ing effect" . They also argued that the QCP focusing 
effect may result in an observable anomaly in the fit de- 
pendence of p/p ratio [2lj . which can be employed as a 
robust signature of the QCP. However, whether the focus- 
ing effect can effectively result in an observable anomaly 
in the fit dependence of p/p ratio has not been worked 
out yet; neither has the corresponding mechanism. In 
this letter, we will formulate the QCP focusing effect on 
the f3t dependence of p/p ratio and apply the UrQMD 
and THERMINATOR model to calculate the dependence 
patterns for three gedanken focused isentropic trajecto- 
ries with different focusing degrees on the QCD phase 
diagram. 



The UrQMD model [22J used here is based on the 
quark, di-quark, string and hadronic degrees of freedom 
and relativistic Boltzmann transport dynamics. It in- 
cludes 50 different baryon species (nucleon, hyperon and 
their resonances up to 2.11 GeV) and 25 different me- 
son species. It is usually used to describe the freeze-out 
and breakup of the fireball produced in relativistic heavy- 
ion collisions into hadrons. The model has successfully 
been applied to reproduce the experimental results from 
SIS/GSI to SPS/CERN energies (H. The THERMINA- 
TOR model [2J| is a Monte-Carlo event generator de- 
signed for studying of particle production in relativistic 
heavy ion collisions from SPS to LHC energies. It imple- 
ments thermal models of particle production with single 
freeze out. The input parameters are those thermody- 
namical parameters at freeze out, such as temperature 
T, baryon chemical potential ps, etc.. 
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FIG. 1: Three gedanken isentropic trajectories on the QCD 
phase diagram with different focusing degrees. All trajecto- 
ries meet at the same QCD critical point (550,159) MeV and 
chemical freeze out point (400,145) MeV on the phase dia- 
gram. 



II. FORMULATION OF THE QCD CRITICAL 
POINT FOCUSING EFFECT 

To quantitatively describe the QCP focusing effect, 
three gedanken isentropic trajectories with different fo- 
cusing degrees, passing through the common QCP and 
chemical freeze out point are constructed on the (/zs,T) 
plane and respectively labeled as Trajectory 1, 2 and 3 
in Fig.l. The location of the QCP (fj, c , TJ=(550,149) 
MeV chosen here is the same as the Ref. [21(, and the 
chemical freeze out point (fi c h, T , c ^,)=(440,145) MeV is 
the statistical model fit result to the data of Pb+Pb 40A 
GeV fixed target reactions at SPS/CERN [Hj]. Besides 
the start and end points, three other (yitB,T) points are 
also sampled for every trajectory. As the actual evolu- 
tion time scale of the isentropic trajectory on the QCD 
phase diagram is unknown, the normalized time is used 
and defined by the normalized path length, t — L/L to t- 
For each isentropic trajectory, the L represents the path 
length along the trajectory from the considered point to 
the QCP and L tot is the length along the trajectory from 
the chemical freeze out point to QCP. The colliding sys- 
tem is evolving from the QCP along the isentropic tra- 
jectory to the chemical freeze out point and assumed to 
be thermodynamical equilibrium. Then, the fj,B and T 
of every sampling point are used as the input thermody- 
namical parameters for THERMINATOR model to cal- 
culate the corresponding antiproton and proton numbers 
along the corresponding isentropic trajectory on QCD 
phase diagram. 

The normalized time dependence of the antiproton 
and proton numbers at the sampled points for the three 
gedanken isentropic trajectories are respectively shown 
in Fig. 2, in which the dash lines are the corresponding 
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FIG. 2: The normalized time evolution of the proton(upper 
panel) and antiproton (lower panel) numbers along the three 
gedanken isentropic trajectories. 



3rd-order polynomial fitting results. It is found in Fig. 2 
that QCP focusing effect results in a sharp decrease of 
proton number and a non-monotonous increase of an- 
tiproton number along the focused isentropic trajectory. 
We can also see large discrepancy among the antiproton 
numbers of the three trajectories, whereas there is little 
discrepancy for the proton number. This means the an- 
tiproton is much more sensitive to the focusing degree 
than proton. In Fig. 3, the time evolution of p/p ra- 
tio along the three gedanken isentropic trajectories and 
the corresponding 3rd-order polynomial fitting dash lines 
are shown. The antiproton to proton number ratio (p/p) 
increases monotonously from QCP (t = 0) to chemical 
freeze out point (t — 1), which can be simply explained 
as the decreasing value of hb/T along each focused isen- 
tropic trajectory [2lT |. 

Although we have obtained three gedanken focused 
isentropic trajectories, the quantitative description of the 
QCP focusing effect on (3 t dependence of p/p ratio and 
also the corresponding mechanism are still ambiguous. 
In qualitative analysis, particles with large fit would be 
emitted earlier from the fireball created in HICs, since the 
mean free path generally grows with increasing hadron 
momentum and becomes comparable with the fireball 
size [2l[. In Fig. 4, the transverse velocity ( (3t ) de- 
pendences of the average emission time ( < t em i SS i on > 
) for antiproton and proton in Pb+Pb 40 AGeV fixed 
target reactions are calculated by UrQMD. It is implied 
that the antiproton and proton, which are dynamically 
emitted from the fireball, would show strong fit — t anti- 
correlation during the cooling down of the colliding sys- 
tem. The general /it — t anti-correlation pattern cannot 
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FIG. 3: The normalized time evolution of the antiproton to 
proton ratio for three gedanken isentropic trajectories. 



be reproduced by hydrodynamic inspired model for their 
particular particle freeze out mechanism. The significant 
increase of p/p ratio with the normalized time t shown 
in Fig. 3, together with the strong fi t -t anti-correlation 
indicates the experimental observable, fit dependence of 
p/p ratio, should be enhanced in the lower fi t region and 
suppressed in the higher fit region [2lj . 

The formulation of the QCP focusing effect on fit de- 
pendence of p/p ratio is based on two simple assumptions: 
one is that the numbers of antiproton and proton evolve 
independently along the isentropic trajectories and the 
other is that the numbers of antiproton and proton emit- 
ted at time t are proportional to their corresponding total 
number (p(t), p(t)). The two assumptions can be formu- 
lated as: 



p(t) = h x / N p (fi u t)dfi t 



p(t) = k 2 x / N p (fi u t)dfi t 



(1) 



(2) 



, where the Np(/3t,t) and N p (f3 t ,t) are the two- 
dimensional fit — t distributions of the antiproton and 
proton, respectively. The two fit — t distributions indi- 
cate antiproton and proton are emitted from the colliding 
system with a finite probability for certain emission time 
t and transverse velocity fit- The k\ and k 2 are the two 
unknown constant coefficients. In addition to equ.(l) and 
(2), we also have the boundary condition: 



p(l) = h x / N p {fi u l)Afi t 



p(l) = fca x / N p {fi u l)Afit 



N p {fi tl t)dfi t dt (3) 



N p {fit,t)dfi t dt (4) 



FIG. 4: The average emission time < t em i sa i 0n > dependence 
of the transverse velocity fit calculated by UrQMD for Pb+Pb 
40 AGeV fixed target reactions at mid-rapidity. 



, which means the numbers of antiproton and proton at 
chemical freeze out point (t = 1) are equal to their cor- 
responding sum of emitted numbers along the isentropic 
trajectory. 

The two constants k\ and k 2 can be determined by 
performing integral over the normalized time t on both 
sides of equ. (1) and (2), respectively as: 



fci 



IoP(t)dt _ flmat 



^ f N p (fi u t)dfi t dt 

JoP(t)dt _ JpP(t)dt 
^ f N p (fi t ,t)dfitdt 



(5) 



(6) 



We introduce D p (t) and D p (t) as the numbers of antipro- 
ton and proton emitted at time t along the isentropic tra- 
jectory, respectively. Then, with equ. (1), (2), (5) and 
(6), we have: 



Dp(t)= I N p {fi tl t)dfi t = '^- 



D P (t)= J N p (fit,t)dfi t = ^ 



p(t) _ p(t) 



sip(t)dt 
m 

foP(t)dt 



xp(l) (7) 



xp(l) (8) 



It is found that the D p (t) and D p (t) are only determined 
by the normalized time t dependence of antiproton and 
proton numbers along the isentropic trajectory, respec- 
tively, which are calculated by THERMINATOR model 
in Fig. 2. 

The experimental observable, fit dependence oip/p ra- 
tio, can be calculated as: 



P(Pt) _ fo N p (fi t ,t)dt 
P(A) ^N p {fi u t)dt 



(9) 
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FIG. 5: fit dependence of p/p ratio calculated by UrQMD for 
Pb+Pb 40 AGeV fixed target reactions at mid-rapidity. 



FIG. 6: fit dependence of p/p ratio for three gedanken isen- 
tropic trajectories. The solid lines are used to guide eyes. 



, which depends strongly on the fit — t distributions of 
antiproton and proton. Although, the two fit — t distribu- 
tions are not exactly known, the anti-correlation between 
fit and t is well known. 

To quantitatively analyze and set the benchmark for 
the fa dependence of p/p ratio with the QCP focusing 
effect, we extract the fit — t distributions of antiproton 
{N^{(i u t)) and proton (N^((3 t ,t)) from the Pb+Pb 40 
AGeV fixed target reactions implemented by UrQMD. 
The emission time t here has been normalized. With the 
two fit — t distributions, the corresponding fi t dependence 
of p/p ratio is calculated by equ.(9) and shown in Fig. 5. 
As QGP phase transition and QCP have not been im- 
plemented in UrQMD model, the monotonous increase 
pattern with (i t up to 0.7 in Fig. 5 is thought to be 
normal and without suffering from QCP focusing effect. 

In order to obtain the anomalous fit dependence oip/p 
ratio resulted from the QCP focusing effects, we modify 
the fit — t distributions of antiproton (Np (fj t ,t)) and pro- 
ton ( (fi t ,t)), which has been calculated by UrQMD 
model. The t distributions in the N£({3 t , t) and Np ((3 t , t) 
are replaced by the distributions Dp(t) and D p (t) de- 
rived from equ.(7) and (8), respectively. For the three 
gedanken focused isentropic trajectories in Fig. I, the 
resulted fit — t distributions for antiproton and proton, 
which have been introduced in QCP focusing effect, can 
be calculated as: 

N ? CP ^=Tmtk* DM <10) 

N?°r(M = j0^ t «nM (ID 

Consequently, the QCP focusing effect has been im- 
plemented in the fit — t distributions for both antiproton 



and proton through the equ.(10) and (11), respectively. 
Thus, the fit dependence of p/p ratio with QCP focusing 
effect can be also calculated by equ.(9) with the modified 
fit - t distributions ( A^ CP (/3 t , t), N^ CP {fi t ,t) ). The 
results for the three focused gedanken isentropic trajec- 
tories on the QCD phase diagram are illustrated in Fig. 
6. Comparing the results of UrQMD calculation in Fig. 
5 and the results with QCP focusing effect in Fig. 6, we 
find that they demonstrate completely opposite depen- 
dence patterns, which indicates the QCP focusing effect 
could actually result in an observable anomaly in the fit 
dependence of p/p ratio. In Fig. 6, it is also noticed that 
the higher focused degree of isentropic trajectory leads 
to the steeper fi t dependence of p/p ratio. Obviously, the 
anomalous fi t dependence of p/p ratio can be employed 
as a sensitive QCP signature. 



III. SUMMARY AND DISCUSSION 

We have formulated the QCP focusing effect on fit de- 
pendence of p/p ratio with some reasonable assumptions. 
The quantitative calculations for the three gedanken 
focused isentropic trajectories have been made with 
UrQMD and THERMINATOR models. In fact, the real 
fit—t correlation pattern for antiproton and proton of 
a colliding system in HICs cannot be measured exper- 
imentally, which directly determines the fit dependence 
of p/p ratio. Therefore, the UrQMD model is used to 
provide fit—t correlation pattern and the THERMI- 
NATOR model is applied to calculate the normalized 
time dependence of the antiproton and proton numbers 
along the isentropic trajectories. Then, the anomalous 
fit dependence of p/p ratios have been obtained for three 
gedanken focused isentropic trajectories, which means 
the QCP focusing effect can efficiently result in an ob- 
servable anomaly in the fi t dependence of p/p ratio. We 
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argue that when the isentropic trajectory of the colliding 
system on the QCD phase diagram is passing through 
the vicinity of the QCP and deformed by QCP focusing 
effect, an observable anomaly in (3 t dependence of p/p ra- 
tio, enhanced in low /3 t and suppressed in high (3 t region, 
will be observed. This anomaly may also be reflected 
in the px spectra of antiproton and/or proton. The 
existing antiproton px spectrum for Pb+Pb 40 AGeV 
fixed target collisions measured by NA49 collaboration at 
SPS/CERN exhibits a steeper exponential slope [2ll.[2^|. 
Finally, we propose that it may be helpful to extract the 
anomalous structures from the px spectrum of antiproton 



and/or proton by performing inverse Laplace transform 
on the spectrum [271 ] . The excitation function of struc- 
ture variable, which should be predefined, can be useful 
to search for the QCP. 
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